Neutron-Proton-Conversion Acceleration at Subphotospheres of Relativistic Outflows 
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We study a type of particle acceleration that operates via neutron-proton conversion (NPC) in in- 
elastic nuclear collisions. This mechanism can be expected for relativistic shocks at subphotospheres 
if relativistic outflows contain neutrons. Using a test-particle approximation, we numerically calcu- 
late the energy spectrum and the efficiency of accelerated particles, and show that a good energy 
fraction of the nucleons may be accelerated. The NPC acceleration may especially be relevant if the 
shock is radiation-mediated, and it would enhance the detectability of GeV-TeV neutrinos. 



PACS numbers: 98.70.Sa, 95.85.Ry 98.70.Rz 

Cosmological gamma-ray bursts (GRBs) accompany 
relativistic jets, where luminous prompt 7 rays are gen- 
erated. Though its emission mechanism has been a long- 
standing problem (e.g., [1]), internal shocks have been 
thought to play an important role [2|. Given that the jets 
are inhomogeneous and have a variety of Lorentz factor, 
a rapid shell finally catches up a slow one and collides. 
In particular, if such collisions happen above the photo- 
sphere, at which the Thompson optical depth tt ~ 1, 
the prompt 7 rays may be explained by optically-thin 
synchrotron emission from electrons accelerated at col- 
lisionless shocks |2|, although this classical scenario has 
several difficulties e.g., the radiation efficiency and the 
inconsistency with the observed spectra [3|, |4| . 

The photospheric scenario is one of the plausible al- 
ternative possibilities, where the prompt emissions are 
produced around the photosphere [a, l6|. In this sce- 
nario, the prompt 7 rays are attributed to modified ther- 
mal [5|-l9| and/or synchrotron emissions [lOllllI- This sce- 
nario could overcome shortages of the classical internal- 
shock model. An interesting channel of the subphoto- 
spheric dissipation exists in neutron-loaded outflows |12| - 
14j . where the hadronuclear reaction between protons 
and neutrons plays an important role and resulting cas- 
cades with Coulomb heating may help to form observed 
spectra [8|. Such neutron-loading is a natural conse- 
quence given that the jet is launched from an extremely 
dense, hot region where the electron capture proceeds [15 1 
(and see [l| for other viable models.) 

In the coming years, neutrino astronomy may pro- 
vide a breakthrough. Assuming that protons are also 
accelerated at internal shocks, PeV neutrinos were pre- 
dicted in the classical scenario [la], and IceCube 17 1 
has constrained reasonable parameter ranges from the 
non-detection of TeV-PeV neutrino counterpart 18l-l21|. 
Different predictions for the photospheric scenario were 
also made Uj, l22|, |23|. Without invoking non-thermal 



sion of IceCube, DeepCore [25|, is able to detect such 
quasi-thermal neutrinos |26l . |27| . Such subphotospheric 
neutrinos are a key to probe physics deep inside the pro- 
genitor. 

Since effective areas of IceCube and DeepCore be- 
come significantly small at lower energies, high-energy 
neutrinos are very important in terms of detectabil- 
ity. However, deeply below the photosphere (tt S> 1), 
there seems a theoretical difficulty in proton accelera- 
tion and the resultant high-energy neutrino production. 
Even if internal shocks can occur, the conventional Fermi 
acceleration would be inefficient at radiation-mediated 
shocks 28|, |29|. The neutron-proton conversion (NPC) 



protons, the inelastic collision model naturally predicts 
multi-GeV neutrinos [12|, |2J|, and the low-energy exten- 



acceleration, which we here explore, can operate even in 
such situations with a neutron-loading. Using Monte- 
Carlo simulations, we demonstrate that a reasonable en- 
ergy fraction of the neutron-flow energy is converted to 
non-thermal nucleons by this NPC acceleration. Here- 
after we use Qx — Q/IO^ in CGS units. 

NPC acceleration: a slow slugger. — The NPC acceler- 
ation was proposed by Derishev et al [30|. In the course 
of conventional Fermi shock acceleration |3ll435[. where 
charged particles cross the shock from the upstream to 
the downstream and come back after being defected by 
magnetic flelds, they implemented conversions between 
charged and neutral particles, e.g., NPCs. We for the 
flrst time perform numerical calculations of this process. 

The advantage of invoking neutrons is highli ghte d by 
considering the possible energy gain per cycle [30|. As 
in the the Fermi acceleration, particles which cross the 
shock experience a Lorentz boost; 

7d ^- 7u = ri.ol7d(l - /Srol/^dMd), (1) 

for the downstream to the upstream, and 

7u ^ 7d = ri.cl7u(l + ^rol/^uMu), (2) 

for the upstream to the downstream. Here, quantities 
subscribed u(d) are deflned in the up(down)stream rest 



frame. We use Fj-ei for the relative Lorentz factor between 
the shock upstream and the downstream, and fj, for the 
pitch angle relative to the shock normal when crossing 
the shock. In a conversion process, a nucleon loses its 
energy as 



7^ «^pn(7- 1) + 1, 



(3) 



with the inelasticity Hp^. As a result, the energy gain per 
cycle can be described as 

(Ef/E,) « Kpn'^^-r.el'Cl - (a*u))(1 + (/id)), (4) 

in the relativistic limit (7^1, Fici ^ 1). Here A'co 
is the number of inelastic collision in the cycle, and ( ) 
means the flux ensemble of the particles who cross the 
shock. One can see that (-Ef/E'i) ex Fioi^ if the particles 
cross the shock after being isotropized, i.e., {fXu} 7^ 1 
and (fid) 76 — 1. 

In the relativistic shock acceleration without conver- 
sions, {Ef/Ei) « 2rroi^ can be realized only in the first 
cycle, and {E{/Ei) < 2 in the successive ones [33l l34j . 
This is because accelerated protons which cross the shock 
from the downstream to the upstream are captured by 
the shock typically in the very early phase of the gyra- 
tion with (//u) ~ 1 — l/Froi^. In the NFC acceleration, on 
the other hand, neutrons that cross the shock can go to 
the far downstream before being converted to protons. 
Then, the converted protons are isotropized in the up- 
stream before being captured by the shock, as long as 
the gyration frequency is much larger than the conver- 
sion frequency (which would be realized in GRB jets) . In 
this case, the NFC acceleration provides a larger energy 
boost per cycle, (Ef/Ei) « 2Kpn^°°Frei^, especially for a 
larger Frei- As a trade-off, the duration of the cycle is es- 
sentially determined by the conversion timescale, which 
makes the NFC acceleration a relatively slow process. 

NPC acceleration in neutron- loaded outflows. — We in- 
vestigate the NFC acceleration in a neutron-loaded out- 
flows. Let us consider an internal shock below the pho- 
tosphere, tt > 1, assuming that two flows contain both 
protons and neutrons. Froton flows are assumed to effec- 
tively form the shock jump in the sense that length scales 
of collisionless dissipation or radiation precursor are rel- 
atively short. Relativistic neutrons go through the flow, 
and they cause inelastic collisions with nucleons with a 
larger length scale 36[ . The conversion channel is domi- 
nated by inelastic nuclear collisions, p + p ^- n + p + TVtt 
and n+p -^ p+p + Nir (see [30| for other cases). For our 
purpose, we can simply assume that (i) a NFC converts a 
nucleon into either of a proton or a neutron with 50% per 
each collision, (ii) the collision is isotropic in the center- 
of-mass frame of incident and target nucleons, and (iii) 
the inelasticity and the cross section are independent of 
the energy, Kpn = 0.5 and Cpn = 3 x 10^^^ cm^, respec- 
tively. Those assumptions can be improved although it 
is not necessary for our present purpose. For example. 



the angular distribution of nucleons can be taken into 
account (e.g., [Ulli]). 

A rapid proton flow with F^ and a slow one with 
Fg collide to form an internal shock at r w 2Fs^ro ~ 
2 X 10^^Fs,2''o.7 cm, where Tq = 10^ro,7 cm is the 
launching radius of the jet. One can estimate the 
Lorentz factor of the shocked region as F w -^/FrFg ^ 
600ri.,3.5^/^rs,2^''^, and the relative Lorentz factor be- 
tween the un-shocked slow shell) and the interaction shell 
as F,ei « 0.5(F/Fs + F,/F) ~ 3Fr,3.5rs,2"^ For a (coast- 
ing) slow flow, neutrons are coupled with protons up 
to the decoupling radius where Tdcc ~ "■ufpn?'dcc/Fs, or 
rdcc ~ LisoO'p„/47rmpC^rs^ -- 5x lO^^Liso.siF^I cm. Here 



Liso/47rFs 



is the baryon number density in 



the upstream rest frame, and L[so is the isotropic lumi- 
nosity of the slow flow. When internal shocks happen 
below rdoc, neutrons are injected into the downstream 
with a Lorentz factor ^ Frei. 

For collisions with the neutron-loaded slow flow at 
Tpn ^ 1, nucleons in the rapid flow are dissipated 
due to inelastic collisions, and the hadronuclear reac- 
tion leads to electrons-positrons, 7 rays, and neutrinos. 
The typical neutrino energy is E^, « 0.1FFj.g[mpC^ ~ 
150rj.p[g 5r2.7 GeV in the observer frame, where ^'^^^ is 
the relative Lorentz factor between the un-shocked rapid 
shell and the interaction shell. Such quasi-thermal neu- 
trinos may be detectable by IceCube-|-DeepCore [2a] . For 
collisions at Tpn > 1, neutrons in the slow shell may be 
dissipated before they reach the un-shocked rapid shell. 
Our goal is to show that a fraction of the injected neu- 
trons re-cross the shock from the downstream to the up- 
stream, and are accelerated up by the NFC mechanism. 

Before proceeding, we should note that our setup in- 
cludes situations where the the conventional Fermi shock 
acceleration would be inefficient. For tt ^ 1 (note tt > 
Tpn), the shock is inevitably radiation-mediated [39|,|40j, 
where the shock width or the deceleration length of pro- 
tons incoming from the upstream to the downstream is 
typically much longer than the isotropization length of 
low-energy cosmic rays 28J, l29|. Such protons cannot 
perceive the enough jump in the flow velocity, which is 
crucial for the energy gain. On the other hand, the mean 
free path of elastic and inelastic collisions must be longer 
than the deceleration length, so neutrons are directly in- 
jected into the downstream with 7d^o ~ Frei- 

Now let us consider a possible acceleration cycle af- 
ter the neutron injection. From Eq.(|3]), the energy gain 
per cycle may be maximized by including the proton 
phase both in the upstream and the downstream to be 
isotropized in the magnetic field, considering the smallest 
number of inelastic collisions. The optimal cycle satisfy- 
ing this criterion is 
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(5) 



0.01 



0.001 



0.0001 




0.1 



0.01 



0.001 



0.0001 



1.5 2 2.5 
log(E [GeV]) 




1.5 2 2.5 3 
log(E [GeV]) 



FIG. 1: The energy spectrum of protons in the downstream for Frei ~ 3 (left) and 5 (right). We set 7d.o = Frci, Tpn = 2, and 
^(1) = 10^. The spectra are normaHzed by the total kinetic energy of the neutron injection. 



(hereafter NPC cycle). The NPC cycle starts from con- 
versions of injected neutrons into protons in the down- 
stream ([n,DS]o => [p,DS]). After being isotropizcd in 
the magnetic field, these protons are re-converted to neu- 
trons ([p, DS] => [n, DS]) while they are advected. A frac- 
tion of the neutrons can cross the shock to the upstream 
([n, DS] — > [n, US]), and again can be converted to pro- 
tons ([n, US] => [p, US]). These protons in the upstream 
are easily captured by the shock, and return back to the 
downstream ([p, US] — ?► [p, DS]). Note here that once the 
protons become non-thermal, the deceleration within the 
shock width can be neglected since such relativistic pro- 
tons and ambient electrons are collisionless. The energy 
gain per NPC cycle is (Ef/Ei) « O-SFi-d^. 

The return probability, Prot, in the NPC cycle can be 
roughly estimated as below. First, both of the two inelas- 
tic collisions must have conversions, which occur with 
a 1/4 chance. Second, only downstream neutrons with 
fid > /3sh,d//3d can cross the shock to the upstream. Here 
/3sh,d is the shock velocity in the downstream rest frame, 
which becomes ~ 1/3 in the relativistic limit. Finally, the 
fraction of neutrons that experience an inelastic collision 
in the upstream is ~ min[l,rpn]. Note that the fraction 
of protons that leave the upstream is quite small for rela- 
tively ordered magnetic fields that we here consider. The 
above arguments yields 



.Prct ~ /npc X — min[l,Tpn] 



(6) 



Here, /npc is a factor to be determined by numerical 
calculations, including all other uncertainties, e.g., the 
fraction of upstream protons which experience inelastic 
collisions before being captured by the shock. 

We can define the acceleration efficacy of the NPC 
mechanism as the energies of accelerated nucleons over 
that of injected neutrons, which is given by £„ 



(1/2), 



'pnj 



Zm/Ei) X Prct) 



JVc- 



pc 



^pn 



Here TVcyc is 



the cycle number, and the pre-factor corresponds the en- 
ergy loss and the survival fraction at ([n, DS]o =^ [p, DS]). 



As we discuss later, N^yc would be at most a few, consid- 
ering other cooling processes. Accordingly, we take only 
the iVcyc = 1 component, which gives J26| 



-npc 



/n 



F 



rcl 



pc 



96 



niin[l,Tpn^] 



(7) 



Note that definitions in [26J are somewhat different from 
those given here |4J]. 

Monte- Carlo simulations. — Here we perform Monte- 
Carlo simulations of the NPC acceleration to justify and 
clarify the estimates above. 

For demonstration, we assume ordered magnetic fields 
parallel to the shock both in the upstream and the down- 
stream, and the compression ratio is the same as the 
baryon density; B^/B^ = n^/n^ — 4(Frci + 3). Note that 
this is not a critical assumption and magnetic fields are 
relevant to isotropize protons. The downstream temper- 



ature can be estimated as Td ~ ("-u^pC'^r^pj/a)^/'* ^ 
1 Piso,52'/Vo,7-i/2Frei,o.5^/'r2.7-i kcV, whcrc Other 
cooling processes than the inelastic nuclear collision can 
be neglected for a few NPC cycles. Consequently, the 
system is parameterized by Fj-ci, Tpn, and ^(1). Here 



C(7) 



i^e,dt 



co,d 



LO, 



,\itco,\n and Wg = 2TTeB/"fmpC 



is the proton-gyration frequency and tco" — napnC is 
the inelastic-collision frequency. When ^(7) S> 1, pro- 
tons are isotropized before the next inelastic collision. 
In the following calculations, we fix ^(1) = 10^, which 
corresponds to a conservative magnetic-field strength of 
Bn ^ 4 X 102 Liso,5irii.3-'r,,2-' G. 

We inject 10^ neutrons setting the initial Lorentz fac- 
tor and pitch angle as 7d,o = rrci and /id.o = — 1, respec- 
tively, and trace the trajectories until the shock sweeps 
the optical depth Tpn, which corresponds to the dynami- 
cal time of the outflow. 

FigH] shows the energy spectra of protons in the down- 
stream normalized by the neutron injection for a fixed 
optical depth, Tpn = 2. The left and right panel shows 
the case of Fid = 3 and 5, respectively. The various 
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FIG. 2: The acceleration efficiency of the NPC mechanism. 
The total energy of accelerated baryons by a single cycle is 
normalized by that of the neutron injection. We fix 7d,o ~ Trd 
and f (1) = 10^. The circles, triangles, and squares correspond 
to Tpn = 0.1, 1, and 2, respectively. 



bumps correspond to the cycle number A^cyc = 0^3. 
The iVcyc = peak is at « 0.5(7d,o — 1) + 1. We confirm 
that the NPC cycle (Eq.®) gives a dominant contribu- 
tion for the A'cyc > 1 components. As expected, the 
peaks are boosted by « CSFrei^ per cycle. 

We also find that, for the Ncyc = 1 bumps, there is 
a non- negligible contribution from a path [n, DS]o => 
[n,DS] -> [n,US] ^ [p, US] -)■ [p,DS]. Thanks to the 
smaller number of the inelastic collision, the energy gain 
in the above path is, in principle, a factor 2 larger than 
the NPC cycle. This path, however, has no proton phase 
in the downstream, and only neutrons which experience 
a large-angle scattering can cross the shock to the up- 
stream. The fraction becomes smaller for a larger Frci 
where most of the scattered neutrons are still directed to 
the far downstream. 

If we approximate the energy spectrum of accelerated 
particles as a power low, dN/dE oc E~'^, the index can 



be estimated as s « 1 -h ln(l/Prct)/ln((^f/£;i)) [3l|, |34| 



(but see also [4l|, |43|). From Eq.(|l]) and ©, one can 
expect the energy spectrum of the NPC acceleration to 
become harder with a larger Fid, which can be seen in 
FigUl although they are not as hard as argued in [30|. 
Fig. [5] shows the acceleration efficiency, where we 

The circles, tri- 

= 0.1,1, and 2, 

for a fixed 



only take the A^cyc — 1 component. 



angles, and squares correspond to r, 
respectively. One can see that e„ 



pn 

ex Td 



Fioi. Also, e„ 



oc 



Fioi^ for a fixed Tpn especially at a 



larger Fj-ei, which is consistent with Eq. ([7]). The en- 
hancement seen at lower Fi-ei is mainly from the path 
[n,DS]o ^ [n,DS] -^ [n, US] ^ [p, US] -^ [p,DS]. The 
results are roughly consistent with the analytical esti- 
mates with /npc ~ 0.1-1. 

Summary and discussion. — We numerically investi- 
gated the NPC acceleration, which may be efficient at 
internal shocks. It may be relevant even for radiated- 



mediated shocks, where the conventional Fermi shock 
acceleration would be inefficient [28|, |29[. This mecha- 
nism may play a role on generating a non-thermal com- 
ponent at subphotospheres of neutron-loaded, relativistic 
outflows. We showed that ~ Frei^min[l,Tpn^]% of the 
neutron-flow energy may be converted to non-thermal 
nucleons with boosts of > O.SFrei^. 

In this work, we adopted a test-particle approximation 
assuming that the neutron fraction is less than unity, 
where the back reaction on the background shock struc- 
ture is neglected. Once the total energy or pressure of ac- 
celerated nucleons becomes signiflcant compared to that 
of the proton flow (rather than the neutron flow) , inelas- 
tic collisions in the upstream contribute to deceleration 
of the proton flow with the length scale « l/n^apa and 
the results should be affected. 

So far, we only took into account the hadronuclear 
collision. In fact, other energy- loss processes may deter- 
mine the maximum energy obtained by the NPC accel- 
eration. In the case of GRBs, the Bethe-Heitler (BH) 
process p + 7— ?>p + e^+e+ would become crucial for 
sufficiently high-energy protons. For a black-body spec- 
trum, this gives a maximum Lorentz factor of 7d,max ^ 
2r7ioC^/CfcB'7d, where C is the pre-factor taking into ac- 
count the effect of the Wien tail. 

In addition, the NPC mechanism becomes inefficient 
for ^(7u(d)) ^ 1: where the pitch angle of a proton is no 
longer isotropized before the next conversion or crossing 
the shock. Then, it becomes difficult to cross the shock 
from the downstream to the upstream. Also, the typical 
pitch angle in the upstream becomes (^u) ~ 1 — l/Prei 
as in the case of the Fermi acceleration, which makes the 
energy gain per cycle negative, {Ef/Ei) < 1, due to the 
inelasticity of the inelastic collisions. This sets another 
constraint of 7d,max ^ ?(!)• Consequently, the maximum 
Lorentz factor by the NPC acceleration can be described 
as 



7d,. 



eSu 



2mcC 



CfceTd ' (TpnmpC^nu 



(8) 



For instance, substituting F — 600, Fioi == 3, rp„ = 1, 
and ^(1) = 10^, which is a possible parameter set for a 
successful GRB jet [26|, the NPC acceleration can give 
7d,inax ^ 10^ if C ^ 10. The by-product neutrino en- 
ergy can be Ei, « 0.05r7dmpC^ ~ 2F2.77d,2 TeV in the 
observer frame. Such s high-energy tail may be cru- 
cial for the detection of subphotospheric neutrinos from 
GRBs [20]. Also, it might affect the signal of hadronic 7 
rays. We should remark that the NPC acceleration may 
operate in failed GRB jets [4^ and proto-neutron star 
winds [36] buried in the stellar material. 
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